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Abstract 

Making use of the Bogoliubov-de Gennes equation, we study the quasi- 
particle spectrum and the vortex core structure of a single vortex in quasi 
2D s-wave superconductors for small pf£,o , where pf is the Fermi momen- 
tum and £o = ■uf/Ao is the coherence length(7i = 1). In particular we find 
that the number of bound states decreases rapidly for decreasing pf£,o- Also 
for pf£,o ~ 1] the Kramer- Pesch effect stops around T/T c ~ 0.3. 
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I. INTRODUCTION 



There are renewed interest in the vortex structure since the discovery of the high T c 
cuprate superconductors. Perhaps the superconductivity of high T c cuprates is characterized 
as (i-wave superconductivity!'! and is very close to the quantum limit 1!. 

Schopohl and Maki have studied earlier the quasi-particle spectrum around a single 
vortex line in terms of quasi-classical equation!'! and predicted a clear four fold symmetry for 
a <i-wave superconductor^. Then a beautiful STM study of vortices in YBCO monocrystals 
was reported!. A big surprise is first of all the quasi-particle spectrum exhibits a circular 
symmetry. There is no trace of four-fold symmetry. Second, there appears to be only a 
single bound state with energy « |A(0) where A(0) (= 260K) is the superconducting order 
parameter at T = OK. A similar bound state energy is observed earlier by a far-infrared 
magneto-transmission from YBCO film by Karra'i et a/.i. If you recall the analysis of bound 
states around the vortex core by Caroli and co-workers0, there should be thousands of 
bound states. Of course in usual s-wave superconductors we have Pf£o ~ 10 3 ~ 10 4 , where 
Pf is the Fermi momentum and £o — v fJ.A(0) is the coherence length(ft = 1). So perhaps 
the single bound state in YBCO suggestsap F £ « 1. At the first sight this suggestion looked 
outrageous, since this implies Ep « 200 ~ 500K in YBCO. By analyzing the spin gap at 
T = OK observed in YBCO monocrystals by inelastic neutron scattering by Rossat-Mignod 
et alM, we can deduce the chemical potential \i given as 
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fi = -3A5(x - 0.45)K (1) 

where x is the oxygen dopage corresponding to YBa 2 Cu30 6+:r . Then for optimally doped 
YBCO we obtain fi = -190K. 

In Dottinger et a/.i the flux flow resistance of 60 K YBCO measured by Matsuda et a/Jil 
is analyzed and theyidentified the Kramer-Pesch effect characteristic to the superconductor 
in the clean limitBa On the other hand apparently the Kramer-Pesch effect is absent in 
90 K YBCO0 implying again perhaps 90 K YBCO is in the quantum limit!. Here Kramer- 
Pesch effect means that the vortex core shrinks with decreasing temperature due to the 
decrease in the occupied bound states around the vortex line. The core size is expressed as0 

f =^^— (2) 
qi A(T)T C ' 1 ; 

This reduction of the core size results to the nonlinear conduct ivityi'lil. Ichioka et a/ill found 
the Kramer-Pesch effect in <i-wave superconductivity with the help of the semi-classical 
approach developed in Ref. [7|. 

But if Pf£,o ~ 1, this implies the semi-classical approach introduced in Ref. [| and Ref. 
|6] is no longer reliable in studying the vortex in the high T c cuprates. So Morita, Kohmoto 
and MakiE-j studied Bogoliubov-de Gennes equation for a single vortex in c?-wave supercon- 
ductors. Indeed choosing Pf£o = 1.33, they can describe gross features of STM result for 
YBCO. On the other hand they have not attempted the self-consistency due to the numerical 
difficulties. Nevertheless they discovered there is a single bound state for pfC,o — 1-33 in the 
vortex of d-wave superconductors. Further there are low energy (E < 0.1 A) extended states 
with four legs stretched in the four diagonal directions (1, 1.0) and (1,-1,0). Recently these 
extended states are rediscovered by Franz and Tesanovicta in a model somewhat different 
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from the one in Ref. [17]. Also in the model of Franz and Tesanovic they don't have any 
bound state. But since they introduced a large short range Coulomb repulsion, it is very 
likely the bound state is completely knocked off. 

In order to avoid the mathematical complication related to d-wave superconductors and 
to prepare for such a study in ci-wave superconductors, we study in this paper the quantum 
limit of a single vortex in s-wave superconductors. Contrary to d-wave superconductors 
the angular momentum around the vortex is a good quantum number, which simplifies the 
problem. Indeed we can follow the formalism set up in Gygi and SchliiterEl except we are 
now interested in the region of small Pf£,o- Recently a related study by Hayashi et a/B is 
reported. Unfortunately they have used the Pf^q dependent cutoff energy rather than the 
constant cutoff energy. So the comparison of systems with different Pf^o appears to have 
involved some artifact. One of our purposes is to find the result free from artifact. 

In the following we limit ourselves to Pf^o = 1,2 and 4 and study the quasi-particle 
spectrum and the shape of |A(r)| as a function of temperature. Even for Pf^o = 1 we find 
3 bound states at T = OK. Also the Kramer-Pesch effect which is present at T/T c ~ 0.5 
appears to stop at T/T c = 0.3 for pf^o = 1- Finally the shape of |A(r)| appears to depend 
crucially on the Uo(r)vQ(r), where uq(t) and vo(r) are the wave function associated with the 
lowest bound state energy. We also show the local quasi particle density of states which 
exhibits remarkable structures. 



II. BOGOLIUBOV-DE GENNES EQUATION 

The spatial dependence of order parameters of the superconductivity is described by the 
Bogoliubov-de Gennes equation. We consider the case with low magnetic field near if cl so 
that there is a single vortex. Therefore we can ignore the vector potential. In this case the 
Bogoliubov-de Gennes equation becomes 

E n u n (r) , (3a) 
E n v n (r) , (3b) 

where u n (r) and v n (r) are quasi-particle wave functions. 

We take z- axis along to the vortex line. We consider nearly two-dimensional case for 
simplicity, where kinetic term associated to the z direction is negligible. In the following, 
we merely consider two-dimensional case and we take a cylindrical coordinate. Taking the 
gauge as A(r) = | A(r)|e , angular momentum of each eigenstate becomes half odd integer 
m + ^0. Then u n (r) and v n (r) becomes as follows, 




Following Gygi and SchliiterEj, we apply Fourier-Bessel expansion with basis; 
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V 2 - /ij u n (r) + A (r) v n (r) 
V 2 - fi) v n (r) + A* (r) u n (r) 
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a/2 r 

-K-Jm+lVOtjjn) It 

to u nm {r) and ^(r), where aj m is j-th positive zero of the Bessel function of m-th order 
J m (x). Then wave functions become as, 

Unm(r) = ^2u nmj (j) mj (r), (6a) 

3 

Vnm{r) = J2v nmj (j) m+lj (r). (6b) 

(6c) 

Here, the boundary condition is such that wave functions are zero at the edge of the disk 
with radius R. Then the Bogoliubov-de Gennes equation becomes, 



a- ^ 2 



2m e V R 

Oij m+1 



2 

- M 
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v 



nmj 



5^ ^jij^nmji E nm V nm j, (^-0 



2m e V i? 
where A JU - 2 is given as, 

f R 

A iu2 = / 0mii( r ) |A(r)|0 m+ i j2 (r)rdr. (8) 
Jo 

The order parameter is given as, 

| A(r) | = g Y 

(r)<f> m+lh (r)[l-2f(E nm )}, (9) 

\E n \<E c m>0 Jij2 

where g is the interaction constant, E c is the cutoff energy and f(E) is the Fermi distribution 
function. We solve these equations self-consistently. We fix the cutoff energy (E c = 5A for 
the case of largest Fermi momentum) and the interaction constant, and choose the Fermi 
momentum p F so that Pf£,o = 1,2 and 4. 

III. RESULTS 

In the numerical calculation, the radius of boundary is taken as R = 10£o and maximum 
angular momentum and number of zero points of Bessel functions are taken so that all of 
the quasi-particle states within cut-off energy are taken into the calculation. 

A. temperature dependence of order parameter 

The temperature dependence of the order parameter is shown in Fig.[l]. Near the bound- 
ary r/£ = 10, the order parameter goes to zero and there is a peak before that. These are 
effects from the boundary condition and the finite size. But core structures are not affected 
by the boundary condition. 
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For pf£,o — 4 at low temperature (T/T c < 0.1) there is a shoulder in |A(r)| at the vortex 
core. This feature comes from bound states in the vortex core. This can be seen from FigfJ|, 
where contributions of the scattering states and the bound states are shown separately. For 
Vf^o = 1, the peak position of the contribution to the order parameter is located slightly 
outside of vortex core. Therefore the core structure is not so much affected by the bound 
states. For larger Pf£o(~ 16), Hayashi et a/B showed the oscillation of the order parameter 
at the inside of the vortex core. It appears that the origin of this oscillation is the reflection 
of the bound state wave functions contrary to their interpretation. 

Slightly increasing temperature, the bound state contribution decreases rapidly and scat- 
tering state contribution remains almost the same. In this temperature range, the core size 
is dominated by the bound states. Above this temperature region, scattering state contri- 
bution decreases with increasing temperature as the behavior of the order parameter of the 
uniform solution. 

We also plot the quasi particle wave function u n ^{r) and v nfM (r) of three bound states 
for Pf£,o = 4 in Fig.^. From these figures, it can be seen that main contribution to the core 
structure comes from the lowest energy bound state. Also we can see u(r) of the lowest 
energy bound state behaves like s-wave function and v(r) behaves like p-wave function. 
Similar behavior also can be seen for second and third lowest bound states. They belong to 
the angular momentum m and m + 1 state respectively. 



B. Quasi-particle spectrum 

In Fig.|| we give quasi-particle energies with angular momentum m + | and we also plot 
those of the uniform state in order to compare the vortex state and the uniform state. For 
small angular momentum, lowest energy states become bound states lowering its energy and 
slightly increasing energies of higher energy states. For larger angular momentum, the quasi- 
particle energy is almost same for vortex and uniform states and there is no bound state. 
Although it is difficult to say the boundary of these two cases in the angular momentum, the 
number of bound states may be counted as 20 for pf£,o = 4,7 for pf£,o = 2 and 3 for pf^o = 1- 
The temperature dependence of quasi-particle energy is shown in Fig|| In this figure we 
plot quasi-particle energies of the several bound states normalized with the order parameter 
A(T) at r = 5.5£o- From this figure we can see that for smaller Pf£,o and higher energy 
states, the quasi-particle energy varies parallel to the order parameter with temperature. 
But lower energy states for larger pp£o, the quasi-particle energy increases rapidly with 
decreasing temperature and becomes constant at low temperature. This behavior of the 
energy of the lowest bound state was already noted by Gygi and SchliiteiS for much larger 
PfCo but the energy increases with decreasing temperature even at low temperature because 
of the large Pf£,o- 



C. Core radius 



As mentioned in [III A| , the core structure depend on the bound states. As defined by 



Kramer and PeschB we calculate core radius £i as 



5 



i = lim^, (10a) 

= f- E E^iWi-2/(iWW°)^(°)- (iob) 

E nm <E c jij 2 

We plot total £f 1 and a contribution from scattering states to £f 1 in Fig.|| The tempera- 
ture dependence oL£l is almost linear in the intermediate temperature region according to 
Kramer and PeschG-a. For pf£,o = 4, there is a substantial decrease of the core radius at low 
temperature because of empty bound states. Although for smaller pf£,o there is a bound 
state contribution, its peak position of the contribution to the order parameter is nearly at 
the edge of the core. Therefore the core radius is not much affected. But for large Pf£,o, the 
peak of the contribution of the bound states are well inside of the core, so the core radius 
shrinks significantly. 

Comparing with Fig. 5 in Ref. [19], our result for small Pf£,o shows saturations of ^f 1 at 



T/T c « 0.1 for p F £o = 4, T/T c » 0.2 for p F £ = 2 and T/T c « 0.3 for p F £o = 1. This comes 
from finiteness of the bound states. 



D. Current density and magnetic field 

Current density is calculated from 

j(r) = ^— E {f(E nm )< m (r)Vu nm (r) + [1 - f(E nm )\ v nm (r)Vv* nm (r) - h.c.} . (11) 

There is a rotational symmetry, so the current has only 9 component; 

je(r) = —Y,{f( E nm)— \u nm {r)\ 2 - [1 - f(E nm )} m + l \v nm {r)\ 2 \ . (12) 

We show this current density at T = 0.1T C in Fig.|^ for each Pf£,o- The peak position of the 
current density is same as the peak position of the lowest energy bound state. Therefore 
smaller Pf£,o the peak is located further from the vortex core. Therefore there is no simple 
relation between the core size and the peak position of the current. 

From Maxwell equation, we calculate magnetic field parallel to the z-axis as, 

4ir r R 

H z (r) = — j e (r)dr, (13) 

C Jr 

and we normalized it with o /27rA 2 and show in Fig.||, where A = (m e c 2 /Aime 2 ) 1 ^ 2 is the 
penetration depth and 0o = hc/2e is the flux quantum. From this figure, we find that for 
smaller Pf£,o the distribution of magnetic field is more extended to outside of the vortex 
core. 



E. Local density of states 



In order to compare with the STM experiments, we calculate the local density of states 
with thermal average, 



6 



N(r, E) = Y.{\ u nm{r)\ 2 f{E nm - E) + \v nm {r)\ 2 f {E nm + E))}. (14) 

nm 

We also take spatial average with a Gaussian distribution with standard deviation 0.1£o- 
This is shown in Fig.|| In this figure there is a oscillation at E/Aq > 1.0, where Ao is the 
average value of the order parameter at the outside of the vortex core. This feature comes 
from the finite size of our system and it is an artifact. 

In addition to this, there are several peaks due to the discrete bound states for energy 
less than the energy gap {E< Ao) and there is a particle-hole asymmetry in contrast to the 
result of Gygi and Schluter&3. Wangand MacDonald" showed this asymmetry by solving 
lattice model and also Morita et a/E and Hayasi et a/.0 pointed out it previously. For 
smaller Pf£,o, this discrete structure is more apparent. 



IV. CONCLUSION 

In summary, we have solved the Bogoliubov-de Gennes equation for a single vortex state 
in s-wave superconductors for small Pf£,o and obtain the electronic structure around the 
vortex core. We show the effect of the discrete bound states to the order parameter structure 
, the local density of states, the current density and the magnetic field. For smaller pf£,o, 
the contribution of the bound state decreases and depends on the peak positions of wave 
functions of quasi particles. Also the discreteness of the bound states is more apparent for 
smaller pf£,o- 

Recently Sonier et a/.il claimed that they obtained the core size of the vortex in NbSe2 
by measuring the peak position of the current density. But there are a few problems in their 
procedure. First, as we have shown there is no simple relation between the core size and 
the peak in the current density. Second even if such a relation exists the relation between 



these quantities are different in the clean limit from the dirty limit result shown in Ref. |22. 
In particular the core size in the clean limit is much smaller than the one in the dirty limit 
with respect to the peak position in the current density. At least they should have analyzed 
their data in the light of the clean limit calculation by Gygi and SchluteitHl. 

Our calculation has been done on the s-wave superconductivity in the quantum limit. 
For high T c cuprates, ci-wave superconductivity is required. Also in the quantum limit, the 
effects from the constraint of the number conservation may be large as pointed by van der 
Mareli, although we have fixed Pf^o f° r temperature region. Therefore, the anisotropy 
of the superconductivity and the constraint of the number conservation must be taken into 
account in the study of the vortex in high T c cuprates. These are left as future problems. 
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FIGURES 



FIG. 1. The temperature dependence of the order parameter A(r) for pf^o = 4 (a), pfCo = 2 

(b) and pf^o = 1 (c). The order parameter is normalized with Ao = A(r = 5.5£o,^ = 0) for each 
Pf£o- 

FIG. 2. Contributions from the bound states and the scattering states to the order parameter 
for p F £o = 4 (a), p F £ = 2 (b) and p F £o = 1 (c), where A = A(r = 5.5£ , T = 0). 

FIG. 3. The quasi particle wave functions u(r) and v(r) and their product with thermal factor 
for the first (a), the second (b) and the third (c) of the lowest energy bound states. 

FIG. 4. The quasi particle spectrum of the vortex state(O) and the uniform state(dashed line) 
for p F £o = 4 (a), p F ^o = 2 (b) and p F £o = 1 (c). 

FIG. 5. The temperature dependence of the quasi-particle energies of the bound states in the 
vortex state. The energies are normalized by the order parameter A(T) at r/£o = 5.5. (a), (b) and 

(c) are for pfCo = 4, Pf£,o = 2 and pf£,o = F respectively 

FIG. 6. The temperature dependence of ^f 1 for pf£,o = 1,2 and 4. Also contributions from 
scattering states for each pf£,o are plotted. 

FIG. 7. The current density at T = 0.1T C for pf£,o = 4 (a), pfCo = 2 (b) and pfCo = 1 ( c )- 
The current density is normalized with its maximum value. The order parameter is also shown for 
comparison. 

FIG. 8. The spatial dependence of the magnetic field H at T = 0.1T C for pf£,o = 1,2 and 4. 

FIG. 9. The local density of states N(r,E) at T/T c = 0.1 for p F £o = 4 (a), p F £o = 2 (b) and 
Pf£o = 1 ( c )- N(r, E) is normalized with the density of states of the normal state Nq. 
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